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We report on studies of the magnetoresistance in strongly correlated 2D electron system in Si in the critical 
regime, in the close vicinity of the 2D metal-insulator transition. We performed self-consistent comparison of 
our data with solutions of two equation of the cross-over renormalization group (CRG) theory, which describes 
temperature evolutions of the resistivity and interaction parameters for 2D electron system. We found a good 
agreement between the p(T, By) data and the RG theory in a wide range of the in-plane fields, 0-2.1 T. This 
agreement supports the interpretation of the observed 2D MIT as the true quantum phase transition. 



PACS: 73.40, 73.43 

Low-tcmpcrature transport in many high mobility 
2D electron systems was found to manifest a critical 
behavior [1-4]. This phenomenon takes place at low 
carrier densities, therefore the electron-electron correla- 
tions play a crucial role. In transport studies, as tem- 
perature decreases, the overall picture of the tempera- 
ture dependence of resistivity p shows clearly distinct 
behavior in two different domains on the density scale. 
At low densities, n < n c , resistivity exponentially in- 
creases with cooling, whereas at high densities, n > n c , 
resistivity significantly decreases (here n c is a sample- 
dependent critical density, ~ 10 11 cm -2 for high mo- 
bility Si samples). In the former domain, transport is 
consistent with a conventional picture of hopping con- 
duction, that is typical for an insulator [5] . In the latter 
domain, far from the critical density (and at p <C h/e 2 ), 
temperature dependence of conductivity was experi- 
mentally shown [6] to be explained by the Fermi-liquid 
effects that were calculated within the framework both 
of the quantum interaction corrections [7] , and temper- 
ature dependent screening [8]. It was so far unclear 
whether the metallic type transport persists to T = 
and whether the cross-over from the metallic to the in- 
sulating behavior (that is observed at low though finite 
temperatures) signifies a true quantum phase transition. 

The successful comparison with the theory of inter- 
action corrections in the high density regime, n ^> n c , 
encouraged us to extend the comparison to the critical 
regime of lower densities n « n c and higher resistiv- 
ities p ~ h/e 2 . The method that now commonly in 
use for this regime and in the diffusive interaction limit 
Tt <C 1 is a generalization of the nonlinear er-modcl the- 
ory, which has been developed by Finkclstein [9]. The 



RG equations [9,10] describe length scale (temperature) 
evolutions of the resistivity and interaction parameters 
for 2D electron system in the first order in p/(irh/e 2 ) 
and in all orders in interaction. 

Earlier [11], only experimental p(T)-data [12] in zero 
field has been compared with one of the RG-cquations, 
while temperature dependence of the interaction pa- 
rameter 72 (J 1 ) [9] was not tested. In this paper we 
have extracted, for the first time, 72 (I 1 ) from the low- 
field magnetoresistance measurements, using cross-over 
RG (CRG) equations, proposed in Ref. [13] for in-plane 
magnetic fields, varying from vanishingly low to large. 
We have compared 72(2^) with theoretical dependence 
calculated from the RG-theory [11] in zero field. Also, 
we have compared our p(T, B») data with the solutions 
of the CRG equations. For these purposes, we mea- 
sured p(T, i?||) in a wide range of the in-plane fields 
£?jl = — 2.5 T. We have found a good agreement of the 
measured p(T, Bit) and 72 (T) with the RG theory. 

Measurements were performed with a Si-MOS sample 
(peak mobility 3m 2 /Vs at T = 1.3 K) of the rectangu- 
lar geometry 5x0.8 mm 2 . We used four-terminal ac- 
technique at 5 Hz frequency; source-drain current was 
chosen low enough, I = 10 11A, in order to avoid elec- 
tron overheating. The sample was located at a precise 
rotation stage that enabled to align its plane parallel 
to the magnetic field directions with accuracy of ~ V . 
The alignment was controlled by observation of the sup- 
pression of the weak localization (WL) peak in p xx in 
the magnetic field B. For studies we have chosen the 
temperature range (1.3 — 4.2) K, because for the stud- 
ied high-mobility Si-MOS sample in the critical regime 
n ~ n c the resistivity exhibits a well-pronounced max- 
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imum at about T max ~ 2 — 3 K with relatively low re- 
sistivity /9 max ~ h/e 2 . For lower densities, the p{T) 
maximum shifts to lower temperatures and higher resis- 
tivities; as a results, /0 max (e /irh) is no longer a small 
parameter as needed for comparison with the one-loop 
RG-theory. 

Figure [T] shows typical temperature dependences of 
the resistivity in the critical regime for various Sii fields 
in the range from to 2.5 T. For the studied sam- 
ple, the metal-insulator transition in zero field takes 



place at the critical density 



0.86 x 10 



li , 



temperatures (see Fig. [T]). At Bit = 2.3 T the maxi- 
mum completely vanishes; for higher fields resistivity 
monotonically increases with cooling that signals the 
onset of the insulating state. The cross-over from non- 
monotonic to monotonic temperature dependence of the 
magnctorcsistance occurs in fields gp,sB\\ ~ fcT max in 
agreement with the predictions of the CRG theory [13]. 

Figure [2] shows the magnetoconductivity data Act = 
er(i?||) — ct(0), replotted from Fig. [T] for various fixed 
temperatures and for Sii varying from to 1.1 T. For 
low fields, b = gfiBBu/kT <C 1, the data, as expected, 
are linear in B 2 . When Act(-B) is plotted as a function 
of b 2 , the slopes of the linear \Aa(B)\ curves increase as 
temperature decreases (compare the data with the guid- 
ing parallel dashed lines in Fig. [2]). This proves that the 
data are taken in the diffusive interaction regime [7]. 
Indeed, in the ballistic interaction regime, where the 
magnetoconductance data scale as (B 2 /T) [7], the slope 
would behave in the opposite way, i.e. decrease linearly 
with cooling. The diffusive interaction behavior can also 
be seen from the inset to Fig. [2j where the same data 
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Data shown in Fig. [I] were taken for the carrier density Cl - 
n = 1.075 x 10 11 cm -2 that is very close to n c . There are 
several evidences for the presented data belong to the 
critical regime: (i) the proximity of n to n c , (ii) the non- 
monotonic p(T) temperature dependence with a clearly 
pronounced maximum and with a sharp drop in p(T) 
at temperatures below the maximum, and (iii) the high 
value of the low- temperature resistivity p fts h/e 2 , that 
is close to the critical value ~ 2/i/e 2 for this sample. 

Rough estimate of the transport time r based on the 
Drude formula shows that Tr is < 0.3 for all "metallic- 
like" curves (i.e. for the curves with p(T) maximum in 
Fig.[T]), over the whole studied temperature range. This 
indicates that the data shown in Fig. 1 belong to the 
diffusive interaction regime. 

In zero magnetic fields, there is a well-pronounced 
maximum in p(T) at T max ss 3K. Application of the in- 
planc magnetic field gradually drives the system to the 
insulating state. As B\\ field increases, the p{T) maxi- 
mum becomes shallow and shifts progressively to lower 
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Fig. 1. p(T) temperature dependences for various in- 
plane magnetic fields varied in steps of 0. 1 T in the range 
(0 - 2.5) T (from bottom to top) 



for a fixed B\\ = 0.7 T arc plotted versus (1/T) 2 . A 
positive curvature of the data demonstrates that the In- 
dependence is even steeper than 1/T 2 . 

For the diffusive interaction regime, the first order 
interaction corrections (which are applicable for low re- 
sistivity/high density regime) result in the magnetocon- 
ductance as follows [14,15]: 



ct(T, B) = a(T, 0) - cn 2 72 ( 72 + 1) 



(gAggll) 
(kT) 2 



(1) 



where c « 0.029, 72 



-FSI0- 



r J' r 



Fermi- 
liquid interaction constant, and n 2 takes into account 
correct number of triplet terms for the multivalley sys- 
tem [11]. Thus, within the framework of the interaction 
corrections, the slope of the magnetoconductance curves 
depends on the interaction parameter. 

According to Eq. (1), the decrease in slopes of the 
Act(6) curves with temperature in Fig. [2] indicates a 
temperature dependence of the interaction parameter, 
which is one of the goals of our studies. However, Eq. ((TJ) 
is valid only for p <C h/e 2 , i.e. for high density "metal- 
lic" regime, where temperature dependence of the inter- 
action constant should be logarithmically weak [15]. On 
the other hand, our main interest is the critical regime 
of low densities and high resistivities p ~ h/e 2 . 
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Fig. 2. Magnetoconductivity in the critical regime ver- 
sus square of the ratio of the Zeeman to thermal energy. 
Different curves correspond to different temperatures, 
varying, from top to bottom, in the range 4.2 to 1.4 K. 
Dashed parallel lines are guides to the eye. Inset shows 
magnetoconductivity at By = 0.7 T versus 1/T 2 



The two-parameter RG-theory has been developed to 
describe the behavior of the 2D system in the critical 
regime for the zero magnetic field case [9,11]. It pre- 
dicts a strong (power-low) temperature dependence of 
72 (T). This theory was already shown [11] to describe 
qualitatively the nonmonotonic p(T) experimental data 
in the critical regime. In RG theory, the p(T) maxi- 
mum signifies a turnover from the "localizing" to "delo- 
calizing" behavior. The (p m ax, T max ) coordinates of this 
point, therefore, are convenient for selecting the proper 
boundary conditions (integration constants). We have 
calculated numerically the 72 (T) dependence by solving 
the system of RG-cquations [11] in the one-loop approxi- 
mation for the two- valley case [11]. This zero field result 
is plotted in Fig. [3] as a solid curve. 

The experimental observation of such strong 72 (T) 
dependence would be a crucial test of the RG-theory. In 
order to extract the 72(2^) dependence from magneto- 
conductance data, a parallel magnetic field is evidently 
required. If magnetic field is weak in all experimentally 
studied temperature range (gpsBw <C kT) the depen- 
dence of the normalized resistance p(T, B^)/ p(B^) mSLX 
on n v p(B\\) max lnT/T(i?||) max is given by the same uni- 



Fig. 3. Temperature dependence of 72 . Solid line - the- 
oretical curve calculated from the RG-equations [11]. 
The symbols denote the empirical 72 obtained from our 
data using Eq. (fTJ (filled circles), and using Eq. ([2]) 
(empty circles) 



versal curve as that for B11 = [9,13]. Here T(B\\) max 
and p(B\\) max denote the temperature and the resistiv- 
ity maximum values at a given B\\ field. Expanding 
T(B») max and p(B») max to the second order in i?y one 
can obtain a novel expression for the magnctorcsistivity: 



p(T,B\\) _ p(T,0) 
P(B\\) max Pmax 



(kT) 2 



(2) 



Here the function 77(72) depends not only on 72 
but on T(i?||) max , p(-E>|| )max, as well as derivatives 
dT(B\\) max /db 2 and dp(B\\) max /db 2 at b = 0. In no- 
tations of Ref. [13] 



t>(72)=ao(72)e 2F( ^- 2F W ax > 



d In p n 



du ^ 

b (u) 



-F{u)+F{-iZ 



dh- 



d In T m 



db 2 



(3) 



6=0 



To compare our data in the low field limit 6cl with 
Eq. (2), for each temperature we have determined the 
difference Sp = p(T, B\\)/ p(B\\) max - p(T,0)/p max and 
plotted it as function of b 2 . The 17(72) function was cal- 
culated numerically, by solving the CRG equations (17) 
and (18) from Ref. [13] in the b 2 approximation. The 
slope of the resulting Sp(b 2 ) curve enables us to extract 
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the experimental temperature dependence of 72 by us- 
ing Eq. (2). The result is shown in Fig. (3) by empty 
circles. One can see that our data agrees with the the- 
oretical 72 (I 1 ) dependence [11] with no adjustable pa- 
rameters in the wide range of temperatures. We show 
in Fig. 3 the data only down to T = 2.2 K at which the 
universal curve given by the one-loop RG theory starts 
to overestimate p(T, 0)/p max . For lower temperatures, 
it is no longer possible to extract 72 (T) from t> (72) since 
the latter becomes non-monotonic. Presumably, it indi- 
cates the significance of higher-order (in p) terms in the 
RG theory at such low temperatures. 

For comparison, we also present in Fig. (3) the 72 (T) 
dependence (filled circles) determined from our experi- 
mental data by using Eq. (1), i.e. the first-order inter- 
action corrections. There is only a qualitative similar- 
ity between thus determined 72 (T) and the theoretical 
RG-result (solid line). The disagreement is not surpris- 
ing because the first-order (in p) interaction quantum 
corrections are inapplicable for the p ~ h/e 2 case. 

It should be noted that the observed strong growth 
in 72 (T) (see Fig. 3) at first sight is inconsistent 
with almost temperature independent Fq , determined 
from Shubnikov-de Haas (ShdH) measurements [16,17]. 
There are several possible reasons for this apparent in- 
consistency: (i) a finite perpendicular and parallel fields, 
such as needed for the observation of beats in ShdH mea- 
surements, affects strongly the solution of the RG equa- 
tions by reducing significantly 72 (T) , (ii) ShdH measure- 
ments of Fq might be made away of the critical regime 
of electron densities, and (iii) ShdH measurements at 
low densities were made at low temperatures T < 1 K, 
whereas the critical 72(2^) behavior in Fig. 3 is estab- 
lished reliably only for higher temperatures T > 2.2 K. 

In Rcf. [13] cross-over RG equations have been sug- 
gested to describe the transition from weak (b <C 1) 
to strong (b 3> 1) parallel magnetic fields. In order to 
accomplish the comprehensive comparison of our data 
with the theory, we have directly compared the mea- 
sured p(T,Bu) dependences with solutions of the cross- 
over RG equations [13] for various fixed magnetic fields. 
For this purpose, we normalized the p(T) data by its 
maximum value p(-B) max for each magnetic field. The 
comparison is presented in Fig. 4. One can see the p(T) 
data agree with the RG theory not only in zero field 
(as was demonstrated earlier [11]), but also in the wide 
range of B\\ fields. 

In low fields -By < 1.5 T there is a quantitative agree- 
ment between the experiment and the RG theory. In 
fields < 0.7 T, as discussed above, the p{T) data collapse 
onto the universal curve for Bu = with no adjustable 
parameters. For higher magnetic field, the data collapse 
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Fig. 4. Comparison of the measured p(T) data (sym- 
bols) with the theoretical dependences (lines) calcu- 
lated from RG-equations [13] for various fixed magnetic 
fields (indicated next to each curve, in Tesla). Data 
and curves are offset shifted vertically by 0.02, relative 
to each other, starting from the zero field curve. In- 
set shows example of the comparison between the p(T) 
data for Bu = 2T and the CRG equations [13] with 
A = 1/27T (dashed curve) and A = 0. 71/77 (solid curve) 



deteriorates because the system enters the non-universal 
cross-over regime. In theory [13], magnetic field enters 
the cross-over RG equations as 



B = A[l + 72 (T max )] 



hT 



(4) 



where the numerical prefactor A = Xjlix. 

In the cross-over regime for B\\ > 1.5 T, there is only 
qualitative agreement between the theory [13] and the 
data as shown in the inset to Fig. 4. However, by chang- 
ing the numerical factor A, we could bring the theoret- 
ical curves into a good agreement with the data. For 
example, for B\\ = 2 T the calculated curve agrees with 
the data when A equals 0.71/7T (see Fig. 4). In even 
higher fields B > 2.3 T, the p(T) maximum vanishes and 
comparison with the theory is no longer possible. We 
mention that in high fields the comparison between the 
data and the RG theory is limited to not too low tem- 
peratures due to the divergence of p(T) (in theory) at 
finite temperature in the one-loop CRG equations [13]. 
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Finally, we wish to note that earlier [11] only one of 
the RG-equations, for p(T) in zero field, has been tested 
by comparing with experimental data [12]. In this com- 
parison the temperature dependence of 72 (T) was taken 
solely from theory and other RG equations were not 
tested therefore. Recently [18], an attempt was under- 
taken to test both variables, p{T) and J2{T), where the 
zero field p(T)-data was compared with solution of the 
RG-cquation, and 72 was obtained from the magneto- 
conductance data using Eq. (1). As discussed above, 
Eq. (1) describes interaction quantum corrections and, 
strictly speaking, is inapplicable in the critical regime of 
p ~ h/e 2 . Indeed, Fig. 3 demonstrates that the 72 
dependence, obtained in this way, is only qualitatively 
similar to the exact solution of the RG-equations [11]. 

In contrast, in this paper, wc compared self- 
consistcntly both variables p(T) and 72 (T) with 
solutions of two RG-equations [11] and found a good 
consistency between the data and the theory. Moreover, 
we have found a good agreement of the p{T,B\\) data 
with the solutions of the CRG-equations [13] for various 
B|| fields. 

In conclusion, wc performed a comprehensive compar- 
ison of the transport and magnetotransport critical be- 
havior in a wide temperature range, with the cross-over 
RG-equations which take magnetic field into account. 
Specifically, we compared with the theory temperature 
dependences of both, the interaction constant 72 and 
resistivity p. We have found that (i) the experimental 
72 (T) determined from the low-field magnetoresistancc 
grows rapidly as temperature decreases, in agreement 
with calculated 72 (T) [11], and (ii) the calculated tem- 
perature dependence of the magnetoresistance is in a 
good agreement with the experimental p(T, B\\) data in 
a wide magnetic field range. This agreement strongly 
supports the theoretical interpretation of the observed 
2D MIT as the true quantum phase transition. 
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